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Abstract 

We use quartic oscillators system with two degrees of freedom to 
model Yang-Mills classical mechanics. This simple model explains 
qualitatively many features reported in lattice calculation of (3 + 1) - 
dimensional classical Yang-Mills system. The largest Lyapunov expo- 
nent (A) and the thermalization time were numerically evaluated. We 
also show, in our model, that A scales with 4th root of energy density. 
Here thermalization is due to relaxation phenomena associated with 
the color degrees of freedom. From the physical picture of thermal- 
ization, we speculate that the system with coherent fields (flux tubes) 
formed in relativistic heavy ion collisions can relax by chaos and the 
estimated thermalization time can be smaller than 1 fm/c. 
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Quantum chromodynamics (QCD) is a Yang-Mills (YM) theory, nonlinear 
and hence exhibits chaos . This property is used to understand the con- 
finement as well as deconfinement problems in QCD. Sawidy JIJ] extensively 
studied chaos in a homogeneous SU{2) Yang-Mills system and speculated 
that it might be related to the confinement problem. Recently, Muller et 
al. || studied the classical YM system in (3 + 1) - dimension on a lattice 
and showed that it exhibits deterministic chaos. They related YM chaos to 
the thermalization of gluons in the deconfined state, which may be formed 
in relativistic heavy ion collisions (RHICs). They evaluated the maximum 
Lyapunov exponent, A, and found that it scales with energy per plaquette 
and estimated the thermalization time of QGP formed in RHICs. But their 
analysis involves heavy numerical computations. It was also not clear why 
A linearly depends on energy per plaquette. In the light of their results, it 
is interesting to reanalyse the earlier work on Yang-Mills classical mechanics 
(YMCM), studied by Sawidy 

In this letter we study this problem using a very simple model, quartic 
oscillators (QO), which is related to YMCM of SU(2) under certain limit. 
The Lyapunov exponent (LE) of QO was studied by Joy and Sabir and 
it's statistical mechanics and thermodynamics was studied by us ||. Here 
we find that maximum LE, A, scales with 4th root of energy density, which 
can be seen both by a dimensional argument as well as numerically and ther- 
malization time is estimated. This is similar to the result of Muller et al, 
but obtained with a simpler model, less numerical computations and with 
a clear physical picture of thermalization based on the statistical mechanics 
of chaotic system. The coherent field energy, like flux tubes formed in rela- 
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tivistic heavy ion collisions, may be redistributed among various degrees of 
freedom associated with spatial, color, chromo electric and magnetic fields 
due to chaos in YM theory. 

The YMCM Hamiltonian is a SU (2) YM Hamiltonian density in temperal 
gauge, Aq = 0, and assuming that the vector potentials A a (t) are functions 
of time only It is given by 

HYM = \Y.'£ + \g 2 >< > (!) 

a a,b 

where g is the gauge coupling constant and a = 1,2,3 are color quantum 
numbers. The Hamiltonian, Hym, may be rewritten in the form 

H YM = H F s+Tym = ~(x 2 +y 2 +z 2 ) + ~g 2 (x 2 y 2 +y 2 z 2 +z 2 x 2 )+T YM ,(2) 

where Hps is called a fundamental subsystem (FS) of YMCM and Tym de- 
scribes quasi-rotational freedoms. Let us consider a simpler two dimensional 
model (a = 1,2 in Eq. (JJ)) and also without Tym- The corresponding 
Hamiltonian is 

H 2 = \{x 2 + y 2 ) + \ g 2 x 2 y 2 . (3) 

On redefining variables Xi = g x and X 2 = gy, it reduces to 

H = g 2 H 2 = l -{Xl + Xl) + \x\Xl . (4) 

It is very similar to the quartic oscillator (QO) system with two degrees of 
freedom, which has been studied extensively both classicaly and quantum 
mechanicaly and is given by |4j] 

H Q = \{Xl + XI) + { ±^{Xi + X$) + \X\X\ , (5) 
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where a is a parameter. For a = 1 it reduces to Eq.(^). QO system is highly 
chaotic for a = 1 and becomes less chaotic as a decreases as shown in Ref. 

I- 

When we make a connection between YMCM and QO, we may note 
that the variables Xi have dimensions of energy (E). Hq is a constant of 
motion, say g 2 e, and has dimension of energy density. Hence we normalize 
all variables by 4th root of g 2 e (say, E) to make them dimensionless. Thus 

we get, 

H' = \{ql + ql) + {1 ^{<lt + 4) + , (6) 

where H' = 1, qi = Xi/ E and = ^ and r = Et. Now we can define the 
distance between two trajectories, D(r), in phase space, as 

D g (t) 



x E(te-^) 2 + fe-^) 2 ), (7) 

\ i=i 

where primed and unprimed variables describe two different trajectories. Su- 
perscript G refers to the fact that it is the general definition used in text 
books |J and Ref. M. Note that without normalized variables, D(t) is not 
dimensionally correct. Now the maximum Lyapunov exponent is defined as 

,. 1, ,D(t), \ e 
Ai = lim - In — )-f = , 8 

where 

A E =lim -ln|^§|. (9) 
Here D(t, E) is given by 



$3 to - xi) 2 + A - x'^/eA , (io) 



Xe is equal to Ai for E = 1 and for arbitrary E or g 2 e, Xe scales with E or 
the 4th root of e which has the dimension of energy. This is similar to the 
results of Muller et ai, obtained numerically, whereas here it follows from a 
dimensional arguments. We have also numerically verified Eq. (§) with our 
model; this is tabulated in Table 1. 

The maximum Lyapunov exponent is evaluated using the procedure of Joy 
and Sabir |4[] using the Hamiltonian, Eq. (|5|), with a = 0.99, which is close 
to 1. For a exactly equal to 1, numerical calculations are not reliable and 
large numerical error develops which can be seen from the numerical values 
of e, which is a constant of motion. The Hamilton equations of motion, to 
be solved numerically, are 

X, = X 3 ; X 3 = (a - l)Xf/3 - X X X 2 2 ; 

X 2 = X 4 ■ X 4 = (a - l)X 2 3 /3 - X 2 X\ , (11) 

which follow from the Eq. (||). These equations are to be evolved along with 
the equation for their variations Yi = 5Xi, which are given by 

Yi = F 3 ; ^3 = ((" - l)Xl - Xl)Y x - 2X 1 X 2 Y 2 ) ; 

% = Y 4 ; % = ((a - l)X 2 2 - X\)Y 2 - 2X 1 X 2 Y 1 ) . (12) 

Then Xe is evaluated using the Eq. @ for different E and is tabulated in 
Table 1. It is in good agreement with the scaling relation. 

However, it should be noted that the procedure used by Muller to evaluate 
A is not reliable here. They obtained A from the slope of the log-plot of D{t) 
Vs t. But here the slope depends on the initial conditions one uses. As an 
example a plot is given in Fig. 1, where two trajectories, separated initially 
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by small distance, are evolved using the Eq. ( |TT| ) for a = 0.99 with different 
initial conditions but with same energy (E = 1). 

Generally in gases and liquids, thermalization is due to collisions and the 
ratio of Kolmogorov entropy (KS-entropy) to collison frequency is related to 
the change in thermodynamic entropy S |7J . In our case there are no collisions 
but the nonlinearity of Yang-Mills equations drives the system chaotic or 
ergotic and hence the thermalization. KS-entropy is related to the sum of 
positive Lyapunov exponents and here we have only one positive LE and 
hence Sks — \e- Let t t h is the thermalization time, then Sxstth ~ S or 

tth « S/S K s = S(E)/X E oc -L . (13) 

If we assume that our system is a subsystem and is in equilibrium with a 
larger system with average energy density (e) propotional to the Ath power 
of temperature (T), then tth is inversely propotional to T; this is consistent 
with earlier (3 + 1) - dimensional calculations || |3|. In addition there is a 
logarithmic dependence on T because of S, which depends on energy density 
(S oc log(e)) as discussed in Ref. ||. 

(3 + 1) - dimensional YM system, discussed in Ref. @, ||, is a system 
with large degrees of freedom due to spatial dependence, and the thermal- 
ization there occurs, due to equal energy sharing between various modes by 
nonlinear interactions. In our case, we have Yang-Mills fields at a point and 
energy sharing is due to their nonlinear interactions between only 2 degrees 
of freedom. It should be noted that recent studies of systems with fewer 
degrees of freedom (even two degrees of freedom), which are chaotic systems, 
do exhibits thermalization or equipartitioning of energy and other statistical 
properties as discussed in Ref. and reference there in. Here it is chaos 
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which plays the role of collisions in statistical mechanics. In statistical me- 
chanics, indeed we need large degrees of freedom to have a chaotic or ergotic 
motion. But here, even two degrees of freedom exhibits chaos and hence al- 
most ergotic because of nonlinearity in the system. To understand it further, 
we may rewrite the Hamiltonian, Eq. (^), as 

H^^Ei + ED + ^Bj, (14) 

using the defmtion of electric and magnetic fields in terms of field tensor 

G7 = d»A» a - VA!> + ge abc A^ c , 

a, b, c are color indices which take values 1, 2, 3 and Lorentz indices 
yu, v = 0,1,2,3 with metric (1,-1,-1,-1). e a b c is antisymmetric Levi- 
Civita tensor. In our case a = 1, 2 and we use hedge hog ansatz and hence 
we have electric fields Ei, Ei and magnetic field B3. Now thermalization 
means sharing energy between electric and magnetic fields. If we start with 
whole energy, say, in electric field Ex, after thermalization the energy will be 
equally distributed between E\ , E 2 and B 3 . As we discussed in [|5|, whenever 
a system is almost chaotic equipartition of energy takes place and 
/. dH 2 \ I. 8H 2 \ I dH 2 \ I 8HA 

\ I1 ai7/ = \ l2 ai7/ = \ I1 a^/ = \ l2 a^/' (15) 

which is same as 

El) = (El) = (Bl) = \e , (16) 

Here the angular brackets indicate the time average, which is also equal to 
the phase space average by ergotic theorem. Hence t t h gives the time required 
to equilibrate energy between electric and magnetic fields and their compo- 
nents. This is a very important property of YM fields which, for example 
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in relativistic heavy ion collisions, will thermalize the coherent fields formed 
immediately after the collision. This is due to the relaxation phenomena 
associated with color degrees of freedom rather than momentum relaxation. 
Further more, if we extrapolate the results of Ref. on the quartic oscil- 
lators, with large degrees of freedom, to our system with large color degrees 
of freedom or large number of YMCM systems with few color degrees of 
freedom, the energy distribution in any one component of fields may be ex- 
ponential decay with decay constant inversely proportional to temperature. 
In other words, corresponding field distribution is Gaussian with it's width 
proportional to temperature. This is similar to the case of momentum re- 
laxation, where one gets Gaussian momentum distribution of particles in 
statistical mechanics. It implies that, at low temperature, the probability to 
find any field component with low energy is high and hence the large number 
of gluons are with low energy. That is, low energy gluons macroscopically 
(largely) occupy YM system at low temperature which may be related to 
confinement. But at high temperature, field components with higher energy 
increase, but those of low energy decrease. Finally the probability to find low 
energy and high energy gluons will be of the same order. Now YM system 
may be in deconfined state. 

So far our discussion is based on A which is evaluated using Eq. (JTj]). 
However, in Ref. |§, they used somewhat different definition for D(t). They 
defined 



where Ei are the electric field components, which reduces (in our notations) 




(17) 
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to 

D M (t) = | ~ X'])\/g = 2| ^Y&l/g . (18) 



Xi and 1^ are evaluated from the Eqs. (|11|) and (12). Even though D M (t) is 
different from D G (t), we get the same results. If we take square root instead 



of modulus in Eq. fll8|) we get Xe about half of Xe obtained from Eq. (fiJjj ). 
This is clear from the expression for D M (t) where the exponentially diverging 
term appears as linear term, whereas in D G (t) it is quadratic. 

In our numerical calculations, as can be seen from the Table 1, theoretical 
scaling behaviour of A with E is confirmed. We also point out that the 
evaluation of A by the slope measurement of log plot of D(t) Vs t, where 
D(t) is evaluted from two independent nearby trajectories, is not reliable. 
Where as standard method given in text books gives reliable results. This 
is probably because of the fact that LE concept is at the linear perturbation 
level, Eq. (0), and not with full nonlinear variations, as in Ref. |§ (by 
following two nearby trajectories separately). 

In conclusion, using a simple QO model to represent the homogeneous 
SU (2) YM fields, many qualitative properties obtained by extensive numer- 
ical calculation of (3 + 1) - dimensional YM fields can be reproduced. We 
obtained the linear relationship between maximum Lyapunov exponent and 
the 4th root of energy density, E, numerically as well as from a dimensional 
arguments. The slope of the linear relationship is 0.4 which is of the same 
order as that of Miiller et ai. We also obtained an expression for thermaliza- 
tion time and explained the physical picture of thermalization in our model. 
The thermalization time is of the order of l/(Ai_E) ~ lfm/c for g = 1 and 
the energy density e = 1 GeV/ fm 3 and inversely depends on the 4th power 



of energy density. In the case of relativistic heavy ion collisions, it may 
give an estimate of the order of magnitude of time required to redistribute 
the energy, stored initially in coherent fields, among all degrees of freedom 
associated with spatial, color, electric and magnetic fields, etc. 
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Figure Caption 



Figure 1: Plots of In | D ^ E ^ | as a function of t for a = 0.99, E = 1 with 
initial conditions X3 = .01 (curve 1), 0.1 (curve 2), 0.4 (curve 3), 1.0 (curve 
4), 1.393 (curve 5) respectively with all other Xi = 0. 
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